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The localized eigenstates of the Harper equation exhibit 
universal self-similar fluctuations once the exponentially de- 
caying part of a wave function is factorized out. For a fixed 
quantum state, we show that the whole localized phase is 
characterized by a single strong coupling fixed point of the 
renormalization equations. This fixed point also describes 
the generalized Harper model with next nearest neighbor in- 
teraction below a certain threshold. Above the threshold, the 
fluctuations in the generalized Harper model are described by 
a strange invariant set of the renormalization equations. 
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In the extensively studied Harper equation [pj 

ipi+i+ipi-t + 2Acos[27r(kr + 0)]V>i = Etpi, (1) 

with a equal to the inverse golden mean, most of the at- 
tention has been focussed on the onset of metal-insulator 
transition at A c = 1 where the quantum states and the 
spectrum exhibit self-similarity and are multifractal. ||] 
In this paper, we show that the fluctuations in the lo- 
calized wave functions of the model for A > 1 possess 
the same complexity and richness as the critical states. 
The universal self-similar fluctuations at the band edges 
are determined by the strong coupling fixed point of a 
renormalization operator. We solve for this non-trivial 
fixed point and obtain the universal scaling ratio for the 
fluctuations. In addition, the stability of the fixed point 
is analyzed by linearizing the renormalization transfor- 
mation. In particular, we study the pertubation asso- 
ciated with a generalized Harper equation describing a 
two-dimensional electron gas with next nearest neighbor 
(NNN) interaction in the presence of an inverse golden 
mean magnetic flux (see Eq. 12). QH The universality 
class for the fluctuations turns out to be unaltered as 
long as the NNN coupling is below a certain threshold. 
Above the threshold, the situation is lot more compli- 
cated as the renormalization flow seems to converge on 
an invariant strange set. 

The wave function tpi is written as 



where 7 is the Lyapunov exponent which vanishes in the 
extended (E) and the critical (C) phase. The localized 
(L) phase is characterized by a positive Lyapunov expo- 
nent corresponding to the exponential decay of the wave 
function. It is assumed that the phase 4> m Eq. (1) is 
chosen so that the main peak of the wave function is at 
i = so that is bounded. ||j7) For the Harper equa- 
tion it has been shown analytically that 7 = ln(X) in 
the L phase. Therefore, the function rji describing 
the fluctuations in the exponentially decaying wave func- 
tion satisfies the following tight binding model (TBM) 
for i > §: 

-rVi+i + <%-i + 2Acos[27r(icr + 0)]rfc = Erji. (3) 
A 

We study this TBM using our recently developed decima- 
tion approach Jl(]||| where all sites except those labelled 
by the Fibonacci numbers F n are decimated. At the n th 
decimation level, the TBM is expressed in the form 



fn{i)rj{i + F n+ i) = rj(i + F n ) + e n (i)r,(i). 



(4) 



The additive property of the Fibonacci numbers provides 
exact recursion relations for the decimation functions e n 
and /„: 
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1 + Af n (i) 
f n -i(i + F n )f n (i + F n ) 



l+Af n {i) 

A = e„_i(i + F n ) + f n _x(i + F n )e n (i + F n ). 



(5) 
(6) 



In this approach, the critical phase with self-similar wave 
functions is characterized by a nontrivial asymptotic p- 
cycle (with the length p equal to 3 or 6 for the Harper 
model) for the decimation functions. M On the other 
hand, the asymptotic behavior of e„(0) and / n (0) deter- 
mines the universal scaling ratios 



0= lim \ V (F pn+j )/r,(0)\; j = 0,...,p-l. 



(7) 
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We apply the decimation method to study the fluctu- 
ations of the L phase where Eq. (3) provides the initial 
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conditions for the decimation functions e2(i) and f2(i) in 
addition to the trivial conditions e\ = and f\ = 1. As 
/2~1/A<lin the L phase, the recursion relations sug- 
gest that /„ ~ l/A Fn_1 , i.e. f n — > as n — > oo. This was 
confirmed by the numerical iteration of Eqs. (5) and (6). 
Because f n vanishes asymptotically, the scaling ratio Q 
can be directly obtained as the limit of |e pn +j(0)| as n 
tends to infinity. At the band edges, we find that |e„(0)| 
converges to the fixed point .1726 for all A > 1. There- 
fore, the fluctuations in the localized wave functions are 
universal throughout the L phase and are determined by 
the strong coupling fixed point of the system (see Fig. 

1)- 

Taking the limit A — ► oo, the TBM reduces to 

rji-i + v cos[27r(W + 4>)]r)i = erji (8) 

with v = 2 and e = lim^oo E/X. For the lower (upper) 
band edge, e is equal to —2 (2). At the band center, 
e = and Eq. (8) is identical to the quantum Ising 
model in a quasiperiodic transverse field at the onset of 
long range order. |(| Therefore, the fluctuations at the 
band center in the L phase are described by the conformal 
universality class of the Ising model. |j| |e„(0)|, and thus 
the scaling factor £, exhibits asymptotically a universal 
3-cycle 0.2307, 0.5904, and 0.2683. 

In order to solve for the strong coupling fixed point at 
the band edges (the 3-cycle at the band center could be 
analyzed in a similar way), we first notice that Eq. (5) 
can be simplified into 

e„+i(i) = -e n -i(i + F n )e n (i) (9) 

assuming that /„ = f n -i = 0. Therefore, this sim- 
ple recursion relation describes both the asymptotic be- 
haviour as n — y oo for any A > 1 and the limit A — > oo 
where /„ = for all n > 1. By replacing the discrete 
lattice index i by the continuous renormalized variable 
x = (— a)~ n < ia >, where < > denotes the fractional 
part M, Eq. (9) tranforms into the form 

e n+ i(x) = -e n -i(a 2 x + a)e n (-ax). (10) 

The high n limit can be studied by introducing the renor- 
malization operator 

T[u(x), t(x)] = [t(-ax), -u(~ax - l)t(-ax)]. (11) 

The importance of this operator follows from the fact 
that the recursion (10) can be represented as the opera- 
tion by T on the pair [e n —i{— ex), e n (x)\. Moreover, the 
limiting function e*(x) = lim„^ oc e„(x) corresponds to 
the fixed point [e*(— ax), e*(x)] of the operator T. An 
initial estimate of the fixed point can be obtained by 
applying T subsequently on the pair [e2(— crx), 63(0;)] ob- 
tained from Eq. (8). It turns out that it is easier to 
expand l/e*(x) than e*(x) (they both satisfy the same 
fixed point equation) so we solve for the coefficients of 



l/e*(x) by truncating the series and applying the New- 
ton method to determine the fixed point of T. The power 
series is convergent in the domain \x\ < 1, and we can 
obtain better and better estimates for the principal scal- 
ing ratio £ = |e*(0)| by increasing the order of the power 
series. Including terms upto the order x 23 we observe 
that e*(0) approaches -0.172586410945 ±HT 12 in agree- 
ment with the result obtained by iterating the decimation 
equations. fH| 

The linear stability analysis at the fixed point shows 
that the renormalization operator T has the unstable 
eigenvalues ±a~ 2 and ±cr _1 (c -1 is a double eigenvalue) 
and the marginal eigenvalue — 1. In addition, there is a 
set of stable eigenvalues which are powers of the inverse 
golden mean. It should be noted that our renormaliza- 
tion operator resembles the one of Ostlund and Pandit 
jjj for the study of the critical point of the Harper equa- 
tion. Although, in their case t and u are 2x2 matrices, 
the eigenvalue analysis is similar in their and our case. 
The variation of e or v in Eq. (8) leads to an asymptotic 
escape from the fixed point in the eigendirection associ- 
ated with the unstable eigenvalue a~ 2 . In the same way, 
the variation of <f> can be related to the eigenvalue — a -1 . 
We expect that some of the remaining unstable eigen- 
values and the marginal one are not "physical" because 
they represent variations which are inaccessible by using 
a TBM to define the pair [u, t] (see ref. 0). We could fur- 
ther generalize the analysis by considering the direction 
of a finite decimation function /„ in function space but 
as explained previously, the pertubation in that direction 
is expected to be irrelevant (with eigenvalue 0). 

It is a characteristic feature of the strong coupling fixed 
point that all eigenvalues are powers of the golden mean. 
In the critical case, where the spectrum is singular con- 
tinuous, there is a non-trivial eigenvalue associated with 
change of energy. We do not expect such an eigenvalue 
to appear in the localized phase where the spectrum is 
point like. In fact, the appearence of the eigenvalue a~ 2 
can be traced to the dual system ]j| for which A — > 00 is 
the weak coupling limit. M 

We next study the generalized Harper equation 

{1 + acos[27r((j(i + i) + 4>)}}tp i+1 

+{1 + a cos[2n(a{i ~h + 

+2A cos[27r(cri + <f>)]ipi = Etpi (12) 

describing Bloch electrons on a square lattice with near- 
est neighbor (NN) coupling anisotropy A and NNN cou- 
pling a. This model was recently studied using analytical 
and numerical methods H as well as applying the deci- 
mation scheme. [j5| The model was found to exhibit a fat 
C phase provided a > 1 and A < a. In the fat C phase 
the model exhibited various new universality classes: at 
certain specific values of A and a, the self-similarity of 
the critical wave functions was described by higher order 
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unstable (with respect to a change of parameter values) 
limit cycles of the RG equations while for arbitrary values 
of the parameters, the RG flow converged on an invari- 
ant set. This implied that the fractal characteristics of 
the wave functions were not self-similar. Since the cycle 
lengths were as high as 24, the question of which param- 
eter values exhibited a limit cycle and which converged 
on strange set remained somewhat open. 

We now apply the procedure outlined above to study 
the fluctuations in the L phase of the generalized Harper 
equation (12). We again make use of the explicit formu- 
las of the Lyapunov exponents as obtained from ref. |Q. 
Numerical iteration of the decimation equations shows 
that for < a < 1, the fluctuations in the L phase are 
determined by the same strong coupling fixed point as 
for the Harper equation. For a > 1, the decimation 
functions for the fluctuations in the L phase are found 
to flow away from the the strong coupling fixed point. 
However, in analogy with the Harper case, the asymp- 
totic behavior of the decimation functions appears to be 
independent of the value of A throughout the L phase 
and is described by the same renormalization equations 
(10) and (11). That is, for all values of a, the decimation 
function /„ renormalizes to zero. Thus we can focus on 
the limit A — > oo, where E/X tends to ±2 for the band 
edges, and the generalized TBM for a > 1 reduces to (<f> 
is 1/2 for the negative and for the positive band edge) 

{±1 + acos[27rcr(i — -)]}r)i-x = a[l — cos^ttoi)]^ (13) 

Because of this simple TBM form (containing no param- 
eter such as the energy which would be known to limited 
precision), the recursion relation (9) for the decimation 
function e„ can be iterated upto 35 times thereby study- 
ing systems of the size 14930351. This accuracy is partic- 
ularly crucial in order to observe higher order limit cycles. 
Our detailed numerical study reveals various limit cycles 
at certain values of a. Writing a — l/abs(cos(2itr)), 
the limit cycles are observed for the rational values of 
r. For example, for r = 1, the period p of the limit 
cycle is found to be 3, while p = 12 for the values 
r = 1/3,1/6,1/8,1/14,1/18,2/9. Furthermore, we ob- 
serve p = 24 for r = 1/7, 1/12, 1/16, 3/14, 3/16, 3/23 and 
p = 18 for r = 1/17,1/19,3/17,3/19. Since the RG equa- 
tions cannot be iterated more than 35 times (because of 
memory limitations), we cannot see higher order cycles. 
Based on this study, we conjecture that for all rational 
values of r, the RG flow converges on a limit cycle of 
period which is a multiple of 3. However, the correlation 
between r and p still remains a mystery. 

For arbitrary values of a, the decimation functions 
do not converge on a limit cycle. The plot of e n (0) vs 
£n+i(0) is found to converge on an invariant set which 
resemples an orchid flower (see Fig. 2). We conjecture 
that the invariant set of the renormalization operator is 
a universal strange attractor. The periodic orbits cor- 



responding to the rational values of r are expected to 
be dense on the attractor. However, we conjecture that 
the probability to hit a periodic orbit is still zero. The 
confirmation of these ideas by an explicit solution of the 
renormalization equations remains open. 

It turns out that the results at the strong coupling limit 
of the generalized Harper equation shed some light on the 
C phase of the model, jl^] This is due to the fact that for 
a fixed value of a, the period of a strong coupling limit 
cycle coincides with the period of a similar limit cycle at 
the critical line A = a. Analogously, the existence of a 
strange set in the strong coupling limit strengthens our 
previous conjecture [^) on the existence of a similar set 
in the C phase. 

In summary, our decimation studies show the existence 
of a new strong coupling renormalization fixed point 
which controls the universal fluctuations of the localized 
wave functions in the Harper equation. Unlike the trivial 
fixed point of the weak coupling limit, the strong coupling 
fixed point describes a new non-trivial universality class 
for the Harper equation. In analogy with the critical 
phase, these fluctuations are characterized by a univer- 
sal scaling ratio which is determined by the value of the 
fixed point function at the origin. We are able to find 
a power series expansion of the fixed point and examine 
the stability of the fixed point under the renormalization. 
It turns out that a change in the NN or NNN coupling 
is an irrelevant perturbation as long as the NNN cou- 
pling is below a certain threshold. Above the threshold, 
the system is outside the basin of the attraction of the 
fixed point and the renormalization flow is attracted by 
an invariant strange set of the function space. Our stud- 
ies put the localized and the critical phase on the same 
footing. The localized phase could in fact be viewed as a 
fat critical phase. 

In some recent studies of the Anderson problem, a con- 
troversy persists concerning the multifractal character of 
the non-exponential part of the wave functions, (lj-jl^] 
The major problem in various numerical studies is that 
due to limited precision, it is not clear whether the stud- 
ies really describe the properties of an infinite system. 
The Harper equation and its generalization provide the 
simplest class of one-dimensional models exhibiting lo- 
calization where the Lyapunov exponent is known ana- 
lytically. Therefore, these models contain no parameter 
which limits the precision of iteration of the decimation 
functions so that extremely large systems can be studied. 
Our results for the self-similar fluctuations in localized 
states of quasiperiodic systems strengthen the likelyhood 
of fractal characteristics in random or aperiodic systems 
exhibiting localization. Since localization is generic in a 
variety of systems ]ri| including the models of quantum 
chaos [ pi , the results described here open a new avenue 
in the localization theory. 
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FIG. 1. The universal fluctuations of a band-edge eigen- 
state in the L phase of the Harper equation. Note that unlike 
the critical phase wave function, the L phase fluctuations are 
not symmetric about the secondary peaks. 

FIG. 2. A two-dimensional projection of the attractor of 
the renormalization flow in the strong coupling limit A — > oo. 
For about 4000 different values of a > 1, decimation equations 
were iterated 35 times and the first 10 iterates were ignored 
as transients. 
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Fig 1 Ketoja and Satija 
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